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Abstract. The literature on Dedekind sums is vast. In this expository paper 
we show that there is a common thread to many generalizations of Dedekind 
sums, namely through the study of lattice point enumeration of rational poly¬ 
topes. In particular, there are some natural finite Fourier series which we call 
Fourier-Dedekind sums, and which form the building blocks of the number 
of partitions of an integer from a finite set of positive integers. This problem 
also goes by the name of the ‘coin exchange problem’. Dedekind sums have 
enjoyed a resurgence of interest recently, from such diverse fields as topol¬ 
ogy, number theory, and combinatorial geometry. The Fourier-Dedekind sums 
we study here include as special cases generalized Dedekind sums studied by 
Berndt, Carlitz, Grosswald, Knuth, Rademacher, and Zagier. Our interest 
in these sums stems from the appearance of Dedekind’s and Zagier’s sums in 
lattice point count formulas for polytopes. Using some simple generating func¬ 
tions, we show that generalized Dedekind sums are natural ingredients for such 
formulas. As immediate ‘geometric’ corollaries to our formulas, we obtain and 
generalize reciprocity laws of Dedekind, Zagier, and Gessel. Finally, we prove 
a polynomial-time complexity result for Zagier’s higher-dimensional Dedekind 
sums. 


1. Introduction 

In recent years, Dedekind sums and their various siblings have enjoyed a new 
renaissance. Historically, they appeared in analytic number theory (Dedekind’s rj- 
function |De| l. algebraic number theory (class number formulae |Me| l. topology 
(signature defects of manifolds HZ| ). combinatorial geometry (lattice point enu¬ 
meration |Mo| l. and algorithmic complexity (pseudo random number generators 
E). In this expository paper, we define some broad generalizations of Dedekind 
sums, which are in fact finite Fourier series. We show that they appear naturally in 
the enumeration of lattice points in polytopes, and prove reciprocity laws for them. 

In combinatorial number theory, one is interested in partitions of an integer n 
from a finite set. That is, one writes n as a nonnegative integer linear combination 
of a given finite set of positive integers. We showed in BDR| that the number of 
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such partitions of n from a finite set is a quasipolynomial in n, whose coefficients 
are built up from the following generalization of Dedekind sums. 

Definition 1.1. For oq, ..., a^, n S Z, we define the Fourier-Dedekind sum 
as 

^ ^ yn, 

a„(ai,...,ad;ao) ■= “ (i _ A-i)... (1 - A“<^) ' 

Here the sum is taken over all oo’th roots of unity for which the summand is not 
singular. 

In m, Gessel systematically studied sums of the form 

E ’ 

A“ = l 

where i? is a rational function, and the sum is taken over all a’th roots of unity 
for which R is not singular. He called them ‘generalized Dedekind sums’, since 
his definition includes various generalizations of the Dedekind sum as special cases. 
Hence we study Gessel’s sums where the poles of R are restricted to be roots of 
unity. 

In Section |21 we give a brief history on those generalizations of the classical 
Dedekind sum (due to Rademacher IE], and Zagier [^) which can be written as 
Fourier-Dedekind sums. Our interest in these sums stems from the appearance of 
Dedekind’s and Zagier’s sums in lattice point enumeration formulas for polytopes 
[Mol IpI BVi DRj . Using generating functions, we show in Section 01 that gener¬ 
alized Dedekind sums are natural ingredients for such formulas, which also apply 
to the theory of partition functions. In Section 01 we obtain and generalize reci¬ 
procity laws of Dedekind |De| . Zagier [^, and Gessel ini as ‘geometric’ corollaries 
to our formulas. Finally, in Sectional we prove that Zagier’s higher-dimensional 
Dedekind sums are in fact polynomial-time computable in fixed dimension. For 
Dedekind sums in 2 dimensions, this fact follows easily from their reciprocity law; 
but for higher dimensional Dedekind sums the polynomial-time complexity does 
not seem to follow so easily, and we therefore invoke some recent work of m and 
[D^. 


2. Classical Dedekind sums and generalizations 

According to Riemann’s will, it was his wish that Dedekind should get Rie- 
mann’s unpublished notes and manuscripts |R.G| . Among these was a discussion 
of the important function 

rj{z) = [|(1 - , 

n>l 


which Dedekind took up and eventually published in Riemann’s collected works 

m 


Definition 2.1. Let ((cc)) be the sawtooth function defined by 


w - 
0 


((^)) : 


1 

2 


if X ^ Z 
if X € Z . 
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Here {x} = x — [x] denotes the fractional part of x. For two integers a and 6 , we 
define the Dedekind sum as 



Here the sum is over a complete residue system modulo b. 


Through the study of the transformation properties of 77 under SL 2 (Z), Dede¬ 
kind naturally arrived at s(a, 6 ). The classic introduction to the arithmetic prop¬ 
erties of the Dedekind sum is m\ The most important of these, already proved 
by Dedekind |De| . is the famous reciprocity law: 

Theorem 2.2 (Dedekind). If a and b are relatively prime then 

,,, ,, , 1 1 /a 1 b\ 

sia,b) + sib,a)--- + -i^- + - + -j . 

This reciprocity law is easily seen to be equivalent to the transformation law 
of the ? 7 -function |De| . Due to the periodicity of ((a;)), we can reduce a modulo b 
in the Dedekind sum: s{a,b) = s(a mod b,b). Therefore, Theorem 12 . 21 allows us to 
compute s(a, b) in polynomial time, similar in spirit to the Euclidean algorithm. 

The Dedekind sum s{a, b) has various generalizations, two of which we introduce 
here. The first one is due to Rademacher who generalized sums introduced by 
Meyer |Me| and Dieter m 


Definition 2.3. For o, 6 G Z, a:, y G R, the Dedekind-Rademacher sum is 
defined by 


s{a,b]X,y) 


E 

k mod b ^ ^ 


This sum posesses again a reciprocity law: 


Theorem 2.4 (Rademacher). If a and b are relatively prime and x and y are 
not both integers, then 

s(a, b; X, y) + s{b, a; y, x) = {{x)){{y)) + ^ + ^^ 2 (ay + bx) + ^^ 2 (x)^ . 

Here 

B 2 {x) := {x - [a:])^ - {x - [a:]) -b ^ 
is the periodized second Bernoulli polynomial. 


If X and y are both integers, the Dedekind-Rademacher sum is simply the 
classical Dedekind sum, whose reciprocity law we already stated. As with the 
reciprocity law for the classical Dedekind sum. Theorem 12.41 can be used to compute 
s(a, b; X, y) in polynomial time. 

The second generalization of the Dedekind sum we mention here is due to 
Zagier From topological considerations, he arrived naturally at expressions of 
the following kind: 


Definition 2.5. Let oi,... , 0 ^ be integers relatively prime to ag G N. Define 

the higher-dimeusioual Dedekiud sum as 


s(ao, ai,..., Uf^) . 


(_l)d/2 


ao 


ao — l , 

E TTKCLi 

cot- 

fc=l 


• cot 


nkod 

00 
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This sum vanishes if d is odd. It is not hard to see that this indeed generalizes 
the classical Dedekind sum: the latter can be written in terms of cotangents |RGi , 
which yields 


s(a, b) = 


1 Trfc 1 

L cot —cot—= --s(&;a,l) . 

k mod b 


Again, there exists a reciprocity law for Zagier’s sums: 


Theorem 2.6 (Zagier). If uq, ... ,ad £ N are pairwise relatively prime then 


d 

^ ^ , Uo? . • . 7 Uj, . . . , Qd) = (/>(ao 7 ■ • ■ 7 ^d) • 

j=0 

Here (f is a rational function in ao,...,ad, which can be expressed in terms of 
Hirzebruch L-functions [^ . 


It should be mentioned that a version of the higher-dimensional Dedekind sums 
had already been introduced by Carlitz m 

a\kx Odkd \ 


E 

,kci mod ao 


ao 


■) 


Berndt ^ noticed that these sums are, up to trivial factor, Zagier’s higher-dimensional 
Dedekind sums. 

If we write the higher-dimensional Dedekind sum as a sum over roots of unity, 

I -h 1 A“‘^ -f I 


s(ao;ai,...,ad) = 


E 




A“<i - 1 ’ 


it becomes clear that it suffices to study sums of the form 


1 


-0 (A- -1) 


T E 

A“0 = l/A 


1 


(A“<i - 1) ■ 


Zagier’s Dedekind sum can be expressed as a sum of expressions of this kind. On the 
other hand, we consider special cases of the Dedekind-Rademacher sum, namely, 
for n £ Z, 


(a,5;p0)= ^ 


ka + n 


k mod b 

Knuth discovered that these generalized Dedekind sums describe the statistics 
of pseudo random number generators. In [H], we used the convolution theorem 
for finite Fourier series to show that, if a and b are relatively prime. 


( 2 . 1 ) 




A- 




(1-A“)(1-A) 


{?} 


1 

4b 


Here {x} = x — [x] denotes the fractional part of x. Comparing this with the 
representation we obtained for Zagier’s Dedekind sums motivates the study of the 
Fourier-Dedekind sum 

an{ai,...,ad;ao) = “ ^E^ (i _ a^i) ■.. (1 - A“<^) ’ 

a finite Fourier series in n. Gessel m gave a new reciprocity law for a special case 
of Fourier-Dedekind sums: 
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Theorem 2.7 (Gessel). Let p and q be relatively prime and suppose that 1 < 
n < p + q. Then 


E 

AP = l5^A 


A’^ 


(1-A«)(1-A) 


+ 


E 




n 

2pq 


1 

pq 


+ 1 - 


12 


1 

pq 



It is easy to see that the reciprocity law for classical Dedekind sums (Theorem 
113 is a special case of Gessel’s theorem. We can rephrase the statement of Gessel’s 
theorem in terms of Dedekind-Rademacher sums by means of for p and q 

relatively prime, and 1 < n < p + q, 



n/1 1 lA 1 1/p 1 q\ If—If— 

2pq 2 \p^ q^ pq) ^ 12 \q^ pq^ pj 2\ p ] 2|gJ 

We will now view the Fourier-Dedekind sum from a generating-function point of 
view, which will allow us to obtain and extend geometric proofs of Dedekind’s, 
Zagier’s and Gessel’s reciprocity laws. 


3. A new combinatorial identity for partitions from a finite set 

The form of the Fourier-Dedekind sum 

2 _^ A~" 

(ai,..., a,; ao) = - (1 - A“0 '' • (1 - A-) 

suggests the use of a generating function 

£/ \ ,_ ^ _2_ 

1 — Z°‘° (1 — ■ (1 — 

In fact, let’s expand this generating function into partial fractions: suppose, for 
simplicity, that n > 0, and oq, ..., are pairwise relatively prime. Then we can 
write 


/w= E 


A, 


d-l-1 


-A 


E ^a Bk 

z — X (z — 1)^ 

k=l ^ ’ 


n ^ 

E ^k 

zk 


A“o = l^A k= 

The coefficient A\ for, say, a nontrivial oo’th root of unity A can be derived easily: 

A A-" 


Ax = lim(z - A)/(z) = -- 

z—*X 


ao(I-A“i)--'(l-A“") ■ 
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Hence we obtain the Fourier-Dedekind sums if we consider the constant coefficient 
of / (in the Laurent series about z = 0): 


A A 

(3.1) const(/)= ^ ri + -"+ E ri + E(-l)'^'= 




The coefficients Bk are simply the coefficients of the Laurent series of / about 
z = 1, and are easily computed, by hand or using mathematics software such as 
Maple or Mathematica. It is not hard to see that they are polynomials in n whose 
coefficients are rational functions of the Oq, ..., To simplify notation, dehne 



(3.2) 


fc=i 

On the other hand, we can compute the constant coefficient of / by brute force: 
By expanding 



we can see that const(/) enumerates the ways of writing n as a linear combination 
of the qq, ... ,ad with nonnegative coefficients: 

(3.3) const(/) = # {(/co, ...,kd)€ : kj > 0, koao H-h kdaa = n] 

= P{ao,...,a,j} (^) • 

This defines the partition function with parts in the hnite set A := {oq, ... ,a„}. 
Geometrically, pa(ji) enumerates the integer points in n-dilates of the rational 
polytope 


V := {(xo, ...,Xd)& : Xj > 0, xqUo H-h XdUd = l} 


This geometric interpretation allows us to use the machinery of Ehrhart theory, 
which will be advantageous in the following section. We next give an explicit for¬ 
mula for the famous ‘coin-exchange problem’—that is, the number of ways to form 
n cents from a finite set of coins with given denominations uq, ... ,ad- comparing 
EU with yields our central result BDRI . 

Theorem 3.1. Suppose ao,...,ad are pairwise relatively prime, positive inte¬ 
gers. We recall that the number of partitions of an integer n from the finite set of 
Oi ’s is defined by 


P{aQ,...,ad}{'n) ■= {{ko,-■ ■ ,kd) ■■ fcj > 0, fcoao H- 'tkdad = n). 

Then 


d 



where q{ao ,..., Ud, n) is given by 113.211 . 


^After this paper was submitted, general formulas for these polynomials were discovered in 


BGKl . 
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The first few expressions for q{ao,..., ad, n) are 


q{ao,n) = — 

0-0 

n 

g(oo,ai,n) =-h 

Oofll 

q{ao,ai,a2,n) = 


1 1 

-1- 

oo ai 

V? n 

2ao^ii72 2 


+ 


1 




ao 


aoai ao02 0102 

fll 02 


Oo 


Oi 


9(00,01,02,03,77) = 

/ 1 

H— j- (-1 

4 \ 000102 000103 

n / 1 1 1 


02 0102 

6O0O1O2O3 

1 

+ 


O0O2 OqOi 


1 


O0O2O3 


+ 


4 \OoOi O0O2 O0O3 


+ 


Oo 


Oi 


0102 
02 


+ 


010203 

1 


12 Y 01 O 2 O 3 


1 


Oo 


+ 


O0O2O3 O0O1O3 
Oo , Oo , Oi 


0103 

0-3 

O0O1O2 
, «1 


0203 


Ol 


24 YO1O2 01O3 02O3 O0O2 

02 02 02 03 


00 O 3 02 O 3 


+ - 


03 


OoOl 


1 / 1 


8 Voo 


0003 0103 

1 1 1 

-h-h — 

Ol 02 03 


OoOl O0O2 


+ 


03 

0102 


4. Reciprocity laws 

We will now use Theorem rm to prove and extend some of the reciprocity 
theorems stated earlier. We will make use of two results due to Ehrhart for rational 
polytopes, that is, polytopes whose vertices are rational. Ehrhart ^ initiated the 
study of the number of integer points (“lattice points”) in integer dilates of such 
polytopes: 

Definition 4 . 1 . Let 7 ^ C be a rational polytope, and n a positive integer. 
We denote the number of lattice points in the dilates of the closure of V and its 
interior by 

L{V,n) ■.= #{nVr]Z’^) and L{V°, n) := # {nV° n , 
respectively. 

Ehrhart proved that L{V, n) and L{V°, n) are quasipolynomials in the integer 
variable n, that is, expressions of the form 

Cd{n) n'^ -\ -h ci{n) n + co(n) , 

where co,...,Cd are periodic functions in n. Ehrhart conjectured the following 
fundamental theorem, which establishes an algebraic connection between our two 
lattice-point-count operators. Its original proof is due to Macdonald 
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Theorem 4.2 (Ehrhart-Macdonald reciprocity law). Suppose the rational po¬ 
lytope V is homeomorphic to a d-manifold. Then 

L{r°,-n) = {-lfL{V,n) . 

This enables us to rephrase Theorem Id. II for the quantity 

P{ao,...,aa}{'^) '= # {(^0, ■ ■ ■ , kd) € : kj > 0 , fcoOQ +- h kdUd = n] . 

By Theorem 14.21 we have the following result. 


Corollary 4 . 3 . Suppose ao,... ,ad € N are pairwise relatively prime. Then 
P{ao,...,aa}i^) = q{ao,... ,ad,-n) + '^an{ao,... ,dj,... ,ad]aj) , 


j=o 


where q{ao ,..., Ud, n) is given by 13.211 . 


We note that we could have derived this identity from scratch in a similar way 
as Theorem rm without using Ehrhart-Macdonald reciprocity. 

The reason for switching to 

P{ao,...,a,}in) = 0 

for 0 < n < ao Qd, by the very definition of p'^^^ ad}(’^)- This yields a 

reciprocity law: 


Theorem 4.4. Let ao,..., be pairwise relatively prime integers and 0 < n < 
ao + ■ ■ ■ + ad- Then 

d 

^ ^ fjyi (ao , ■ • ■, dj ,..., ad, aj ) — q(ao-i ■ ■ ■, adj n') , 
j=o 

where q(ao ,..., ad, n) is given by 13.211 . 

For c? = 2, oo = p, ai = g, 02 = 1, this is the statement of Gessel ’s Theorem 12.71 
which, in turn, implies Dedekind’s reciprocity law Theorem 12.21 

To prove Zagier’s Theorem 12.til in the language of Fourier-Dedekind sums, we 
make use of another result of Ehrhart [E] on lattice polytopes, that is, poly¬ 
topes whose vertices have integer coordinates. Recall that the reduced Euler 
characteristic of a polytope V can be defined as 

X(V) 

cr 

where the sum is over all sub-simplices of V. 

Theorem 4.5 (Ehrhart). Let V be a lattice polytope. Then L(V,n) is a poly¬ 
nomial in n whose constant term is x(T’). 

We note that the polytope P corresponding to p^ag,...,ad} ("a) is convex and hence 
has Euler characteristic 1. If we now dilate P only by multiples of oq • • - ad, say 
n = ao ■ ■ ■ Odm, we obtain the dilates of a lattice polytope. Theorem im simplifies 
for these n to 

d 

/^{ao,...,ad} (^0 ’ ■ ' ■ ■ ■ 5 ^d-i ^0 ’ ’ ’ ^ ^ ^0 (^0? • ■ • 5 ■ ■ ■ 5 ) 7 

i-0 
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by the periodicity of the Fourier-Dedekind sums. However, x(7^) = 1, and Theorem 
14.51 yields a result equivalent to Zagier’s reciprocity law for his higher-dimensional 
Dedekind sums, Theorem 

Theorem 4.6. For pairwise relatively prime integers oi,..., ad, 

d 

CTo(ao,Ud; Oj) = 1 - g(ao,..., a^, 0 ) , 

3=0 

where q{ao ,..., aa, n) is given by 13 . 211 . 

5. The computational complexity of Zagier’s 
higher-dimensional Dedekind sums 

In this section we give a proof of the polynomial-time complexity of Zagier’s 
higher-dimensional Dedekind sums, in fixed dimension d. In [HP], there is a nice 
theorem due to Barvinok which guarantees the polynomial-time computability of 
the generating function attached to a rational polyhedron. We will use his theorem 
for a cone. First, we mention that a common way to enumerate lattice points in a 
cone 1C (and in polytopes) is to use the generating function 

/(/C,x):= ^ X-, 

mG/CnZ^i 

where we use the standard multivariate notation x™ := x™’ It is an 

elementary fact that for rational cones these generating functions are always rational 
functions of the variable x. Barvinok’s theorem reads as follows: 

Theorem 5.1 (Barvinok). Let us fix the dimension d. There exists a polyno¬ 
mial-time algorithm, which for a given rational polyhedron 1C C 

/C = {x G : (ci,x) < (3i,i = where Ci G and fii G Q 

computes the generating function 

/(/C,x):= ^ x” 

in the form (a virtual decomposition) 

where Ci G { — 1, l},ai G Z, and bn, ...,bid is a basis of Tfi for each i. The compu¬ 
tational complexity of the algorithm for finding this virtual decomposition is 
where C is the input size of 1C. In particular, the number I of terms in the summand 

IS 

Thus Barvinok’s algorithm finds the coefficients of the rational function /(/C, x) 
in polynomial time. In DR| . on the other hand, the generating function /(/C,x) is 
given in terms of an average over a finite abelian group of a product of d cotangent 
functions, whose arguments are in terms of the coordinates of the vertices which 
generate the cone 1C (these are the extreme points of 1C whose convex hull is 1C). 
This is the main theorem in [DRj and we apply it below to a special lattice cone 
which will give us the Zagier-Dedekind sums we want to study. 
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The following theorem is part of a bigger project on the computability of gen¬ 
eralized Dedekind sums in all dimensions. A slightly different proof is sketched in 


m 


Theorem 5 . 2 . For fixed dimension d, the higher-dimensional Dedekind sums 


(— 1 )'^/^ irkai nkad 

s(ao; oi,..., Od) = - cot-■ • • cot- 


ao 


k=l 


ao 


ao 


are polynomial-time computable. 


Prooe. Let 1C C be the cone generated by the positive real span of the 
vectors 


vi = (l,0,...,0,ai) 

V2 = (0,1,0,..., 0,02) 


Vd 

Vd-i-l 


= (0,...,0,l,ad) 

= (0, ...,0,ao) 


Then the right-hand side of the main theorem of [DRJ is in this case 

oq —1 d 


1 


2‘^+^ao 


En 


1 -f coth — (s -I- iojk) 
ao 


k=i j=o 

When we compute the coefficient of in this meromorphic function of s, we arrive 
at the following higher-dimensional Dedekind sum: 


ao —1 

E 


Trkai 

ao 


• cot 


irkad 

ao 


plus other products of lower-dimensional Zagier-Dedekind sums. By induction on 
the dimension, all of the lower-dimensional Zagier-Dedekind sums are polynomial¬ 
time computable, and since the left-hand side of the main theorem is polynomial¬ 
time computable by Barvinok’s theorem, the above Zagier-Dedekind sums in di¬ 
mension d is now also polynomial-time computable. □ 
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